On recombination in strong laser fields: effect of a slow drift 
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The dynamics of the recombination in ultrastrong atomic fields is studied for one-dimensional 
models by numerical simulations. A nonmonotonic behavior of the bound state final population as 
a function of the laser field amplitude is examined. An important role of a slow drift of an electron 
wave packet is observed. 
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I. INTRODUCTION 

It is well known that during atomic photoionization 
in ultrastrong laser fields a peculiar effect occurs: for a 
given pulse length an interval of field intensities exists 
in which the ionization probability is a decreasing func- 
tion of the field intensity. This is called adiabatic atomic 
stabilization against photoionization 1,2]. The motiva- 
tion of the present work is to look for a similar effect in 
the inverse process, namely atomic recombination in ul- 
trastrong fields. The main difference between these pro- 
cesses is due to different initial conditions. The process 
of strong-field recombination has recently been examined 
by Hu and Collins 3], who however investigated the de- 
pendence of the process on the initial momentum of the 
incoming electron. 

During the recombination one should expect dynam- 
ical effects similar to those observed in the case of ion- 
ization. As it has been shown in the works devoted to 
the atomic stabilization, in laser fields of order of a few 
atomic units most of the electron wave packet moves as 
a whole, performing oscillations in the rhythm of the 
field. Additionally, long time oscillations of the packet 
are possible which are due to an asymmetry of the inter- 
action of its different parts with the binding potential; 
in our earlier papers [J, IS 13 we have called them a slow 
drift. The slow drift has the range equal to that of the 
classical oscillations of a free electron in the laser field 
and it was possible to give an analytic formula for its 
frequency in the case of the binding potential being a 
rectangular well. The fact that a significant part of the 
packet moves without changing of its shape, in particular 
of its width, can be explained in terms of the so-called 
Kramers-Henncberger (KH) well, which is the mean po- 
tential (zero-th term of the Fourier expansion) of the os- 
cillating nucleus in the electron's reference frame. The 
stabilization of the wavepacket can be interpreted as its 
trapping in the eigenstate of the KH well. In the labo- 
ratory frame the KH well oscillates in the rhythm of the 
laser field and may also perform the slow drift, which was 
described by a generalized model of Ref. . For longer 



times a decay of the KH state can be observed, due to an 
influence of higher terms of the above-mentioned Fourier 
expansion, and the corresponding decay rate has been 
evaluated ■ The probability of finding the electron in 
the atomic bound state is thus a complicated function 
of time, which reflects a rapid oscillatory motion, a slow 
drift and a finite lifetime of the KH eigenstate. More- 
over, a necessary condition for the stabilization to occur 
is that the pulse have such a shape that the corresponding 
electron classical trajectory should remain in the neigh- 
borhood of the nucleus. 

In this work we present the results of numerical simu- 
lations of the dynamics of recombination for two model 
one-dimcnsional systems: a rectangular potential well 
and a long-range Coulomb-like potential. We will exam- 
ine the time dependence of the bound state population 
(initially equal to zero) and the dynamics of the wave 
packet. As the measure of efficiency of recombination we 
take the final population of the bound state. In particu- 
lar we study how the details of the dynamics depend on 
the amplitude of the laser field and as a consequence on 
the details of the slow drift. 



II. NUMERICAL APPROACH 

Strong-field recombination has been investigated in 
which an initially free electron modeled by a Gaussian 
packet, intially resting at twice the amplitude of classi- 
cal free oscillations, is carried towards the ion in the laser 
field. The ion has been modeled by two kinds of a binding 
potential. 

The dynamics of recombination has been examined 
by numerically solving the time-dependent Schrodinger 
equation in one spatial dimension. The computations 
have been performed on a grid of 16384 points. The 
space covered by our grid extends from —100 to 100 a.u. 
(the grid step is about 0.0122 a.u.). Calculations have 
been performed for two atomic models: the rectangular 
symmetrical well and a soft core atom potential. The 
rectangular well 
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of the depth equal to Vi = 2.049 a.u. and half-width 
equal to ai — 0.122 a.u. is placed in the middle of the 
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simulation space. The function is the Heaviside step 
function. The same well parameters were earlier used in 
the studies of stabilization Q. We have also used the 
model atom potential proposed by Eberly et al. 
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To simulate the evolution of a one-dimensional wave- 
function we use the Crank-Nicholson scheme, which for 
Schrodinger equation leads to the often used symmetri- 
cal Cayley's formula 9J. The simulations have been per- 
formed in the dipolc approximation in the length gauge, 
i.e. the Hamiltonian of the system reads (in atomic units) 



(3) 



where the laser electric field e{t) — eoO{t)Q{tf — 
t) cos{ujt — tt) is directed along the cc-axis and acts only 
during the time interval tf. The period of the laser is 
T = 2TT/UJ, where lu ^ I a.u. is the laser frequency. 
We have changed the laser field amplitude Eq in the 
range between 1 to 15 a.u. The time step was equal 
to At = T/10^. The convergence of the computer simu- 
lations has been checked on the grid with 65536 spatial 
nodes and alternatively with a ten-times smaller time 
step. No important differences have been visible. 

The initial position of the wavepacket, initially taken 
to be Gaussian \tp\'^ = exp[— (x — xo)^/{2S^)] with the 
variance (5^ = 0.25 a.u., has been prepared very carefully, 
independently for each simulation. The main condition 
is that the electron should be placed over the well after 
the first half of the period, i.e. x{T/2) = 0, x being the 
packet's center of mass. Moreover, we require it to have 
at that time zero velocity to make the recombination 
more probable. This is why we have chosen the packet's 
initial position xq = —2eo/uj'^ and initial velocity 
equal to zero so that the first classical turning point 
of a free oscillation occurs over the well. Note that 
the dependence of the recombination efficiency on the 
initial packet velocity has already been the subject of 
investigations of Ref. • 
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FIG. 1: The population of the ground state in the model 
atomic potential in the recombination process as a function 
of the field intensity eo in the case of the field being switched 
off at exactly ST, solid line; at the time instant of a maximum 
population in the vicinity of ST (see comment in the text), 
dashed line; at the time instant of a maximum population in 
the vicinity of 2r, short-dashed line. 
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FIG. 2: As in Fig. 0but for an only bound state of a well 
potential 



III. RESULTS AND INTERPRETATION 

In the case of the recombination process one observes 
a nonmonotonic dependence of the bound state popu- 
lation on the field intensity. It is well visible in Fig. 
where the ground state population is shown as a function 
of the field intensity. In each case the rectangular pulse 
has been switched off at the time instant closest to three 
optical periods (3T) at which the ground state popula- 
tion exhibits a maximum (except for the solid line which 
has been obtained exactly at 3T). At the time 2T the 
maximum is not yet formed. A wide maximum is seen 
for £o w 2.5 a.u. both for the pulse duration of exactly 



3T and approximately 3T. A similar nonmonotonicity 
is observed in the case of a potential well (see Fig. [SJ. 
The figure suggests an existence of even more maxima 
which are however less pronounced for short times. This 
is confirmed for longer pulses in Fig. |31 (upper plot) in 
which one can see two separated maxima, one of which 
occurs for eo ~ 2.5 a.u. for all pulse lengths, while 
the other one, being wider, shifts slowly for increasing 
pulse lengths. The origin of the peaks is different. The 
left peak can be attributed to a long lifetime t of the 
eigenstate of the Kramers-Henneberger well; according 
to Ref. i T exhibits a maximum for such eg for which 
Ji(£o'y2w/w^) = 0, i.e Eq =2.71 a.u. and 4.96 a.u. (for 
the range of field amplitudes examined in this work) . The 
position of the right hand side peak for different values of 
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the field intensity and of the pulse duration is evidently 
connected with the slow drift of the wave packet. A max- 
imum is observed if at the moment of the field switch-off 
the slow drift of the main part of the packet is at such a 
stage that it is located at the well, which does not usually 
occur after an entire number of periods, and has a veloc- 
ity close to zero. A similar field dependence of the bound 
state population can be seen in the case of the ionization 
(lower plot). Again one can see a peak for Eq w 2.5 a.u. 
for all the pulse durations and a second peak the position 
of which changes. One can also notice a small peak at 
about 5 a.u. (second minimum of Ji and thus of l/r). 
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FIG. 3: The population of the bound state of the well poten- 
tial in the recombination process (upper plot) and ionization 
process (lower plot) as functions of the field intensity £o in 
the case of the field being switched off at exactly lOT, 20T 
and 30r. 

In Figs 21 and [3 we show the evolution of the wave 
packet for different values of the field amplitude, for 
the potential well and for the model atom, respectively. 
For relatively small field amplitudes each field oscillation 
causes tearing the packet into smaller parts, typical of the 
regime of multiphoton ionization. For the field amplitude 
being of order of 2.5 a.u. one can already distinguish a 
main part of the packet, represented by the oscillating 
streak. A slow drift is imposed on the rapid oscillations 
with the field frequency u>. For the potential well the 
frequency of the slow oscillations may be evaluated using 



the formula Q — y 2aiVi/{\/2Tr(T^), where a — IeqIlo^ 
is the width of the wave packet trapped in the KH well 
0. Due to the slow drift, after an entire number of field 
periods, the packet ends up at different distances from 
the potential center, depending on the field amplitude. 
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FIG. 4: The time evolution of the wave packet recombination 
in the well potential for different values of the field intensity. 
The level of blackening is proportional to the squared modules 
of the wave function. 
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FIG. 5: As in Fig. |l]but for a model atomic potential. 

The infiuence of the atomic potential is for our data 
stronger then in the case of the potential well, which is 
manifested by a larger frequency of the slow drift and 
a more effective recombination and population trapping 
(darker line) and less population thrown away as isolated 
small portions. In the case of the well potential, for large 
times the streak for eg = 2.5 a.u. is significantly darker 
than for 2.0 a.u. and 3.0 a.u. (not shown), which cor- 
responds to the first maximum of Fig. |21 (upper plot). 
There is also a correspondence between the level of pop- 
ulation shown in Fig. |31 (upper plot) and the position of 
the streak of Fig. 0] For example for t — lOT and eo = 5 
a.u. the streak is located at the position of the well x = 
while for £q — 7 a.u. its distance from the well is close to 
maximum due to the slow drift; this corresponds to the 
fact that the bound state population in Fig. |31 (upper 
plot) for £0 = 5 a.u. is larger than for Eq = 7 For 
t = 30r the situation is reversed due to the slow drift. 

For larger fields, especially in the case of atoms, the 
streak splits into two ones, which corresponds to popula- 
tion trapping in a superposition of two eigenstates of the 
KH well. One can also see a fiow of the population from 
one streak to the other, which corresponds to the popu- 
lation oscillating between those two states. In the case of 
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FIG. 6: The time evolution of the bound state population in 
the case of a well potential for the field intensity eo = 5 a.u. 



the rectangular well^ for which an analytical form of the 
KH well is known 4], we have verified that for the field 
amplitude eg = 7 a.u. there are indeed two eigenstates, 
the energy difference of which is approximately 0.016 a.u.; 
the corresponding period of oscillations agrees with the 
period of the population flow of Fig. 21 

The presence of the slow drift is reflected in the 
time dependence of the bound state population visible 
in Fig. El The maxima of the comb-like structure 
appear twice in each optical cycle at the time instants 
at which the packet's velocity is zero (turning points). 
The heights of the maxima depend on the current stage 



of the slow drift. The frequency of the exchange of 
the two families of the peaks, appearing respectively 
after an even or odd number of half-cycles, is the 
same, both in the case of ionization and recombination, 
as the frequency of the slow drift, visible also in Fig. ^ 



IV. CONCLUSIONS 

A nonmonotonic behavior of the bound state final 
population as a function of the field amplitude has 
been observed in numerical simulations of strong-field 
recombination for one-dimensional models. The maxima 
may have two origins. First, they appear for such laser 
field intensities and pulse durations for which the packet 
ends up in such a stage of its rapid oscillations and of its 
slow drift that it is located at the nucleus and has zero 
velocity. Second, a maximum, located independently 
of the pulse duration, is due to a maximum lifetime of 
the KH bound state. For such fields that more than 
one KH state exist and are significantly populated, one 
additionally observes a third kind of oscillations of the 
wave packet (in addition to the oscillations with the 
optical frequency and the slow drift), which are due 
to the population oscillating between the KH eigenstates. 
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